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In this work a field theoretical model is constructed to describe the sta- 
tistical mechanics of an arbitrary number of topologically linked polymers in 
the context of the analytical approach of Edrwards. As an application, the 
effects of the topological interactions are studied in the one loop approxima- 
tion. A natural way to include in the treatment also more sophisticated link 
invariants than the Gauss linking number is outlined in the Conclusions. 



I. INTRODUCTION 

Long flexible polymer molecules provide an ideal environment to study the effects of 
topology in statistical mechanics and condensed matter physics, but so far a comprehensive 
theory of topologically linked polymers does not exist. Until now, only the physics of un- 
entangled open chains has been understood to a satisfactory extent 0. In this letter the 
problem is tackled in the context of the so-called analytical approach, consisting in a number 
of different methods based on Edwards pioneering work , in which the fluctuations of each 
polymers are treated as a random walk and described via a path integral in the Wiener 
sense. Analytical approaches of this kind have already been applied to several complex sit- 
uations ||- Jl2| and seem very promising in explaining the physics of topological entangled 
polymers. Despite of that, they are affected by two main limitations. First of all, the Gauss 
linking number, a weak topological invariant, is used to distinguish the topological states of 



the system. Unfortunately, more sophisticated knot invariants like those described in [13 



and [Tj] cannot be applied in practice, since they have no immediate relation to the physical 
conformation of the polymer |J. On the other side, one of the advantages of the analytical 
approach is to provide a field theoretical description of a system of polymers, establishing 
in this way a connection to the theory of critical phenomena in condensed matter and high 
energy physics. However, a mapping from the statistical mechanical problem of polymers 
to field theories has been realized only in a limited number of cases, mainly concerning the 
fluctuations of a single test polymer in a background of static polymers or fixed obstacles. 
Remarkable exceptions occur when two polymers are very close f|, Hl2| , a situation which is 
particularly relevant in the study of DNA molecules |H| , or when the number of polymers 



becomes large [p]] . Despite of these successes, the exact analytical treatment of a system 



of N polymers has been a formidable problem already mentioned in the pioneering work on 
topological entanglement of Brereton and Shah || and remained unsolved until the present 
days. 

To solve it, we construct a Chern-Simons (C-S) [16| based model, which generalizes a 
previous work of the authors in which the two polymer case was discussed |17| . With respect 
to ref. [T7J], the main difficulty is that one needs N ^~^ independent topological numbers 
to distinguish the possible topological states of iV polymers according to the Gauss link 
invariant. As a consequence, one has to find a suitable C-S field theory which decouples 



the actions of the polymers as in [|17]] and simultaneously is able to accommodate the huge 
amount of topological configurations allowed by the system. 

The material contained in this paper will be divided as follows. In the next Section 
a topological Ginzburg-Landau model is constructed which describes the fluctuations of a 
system of N topologically linked polymers. The topological relations are specified using the 
Gauss linking number. Section III contains a perturbative study of the case N = 2. In par- 
ticular, the one loop effective potential is computed proving that the topological interactions 
do not influence the critical behavior of the system in this approximation. Moreover, it is 
shown that the second topological moment can be exactly computed in terms of one loop 
Feynman diagrams. Finally, a detailed discussion of the results has been given in Section 



IV. In that Section a natural method to include higher order topological invariants is out- 
lined. To this purpose, a new C-S field theory is defined, which is abelian but exhibits cubic 
interactions among the C-S fields. In this way such theory can generate radiative corrections 
containing knot invariants with the same mechanism acting in the case of nonabelian string 
theories P|, M. 



II. THE iV-POLYMERS PROBLEM 

Let P 1 ,...,P N be a set of topologically linked polymers, describing in the space the 
trajectories Tj, i — 1, . . . , N defined as follows: 

Ti = {rj(sj) |0 < Si < Li ; r;(0) = r-, r^Lj) = r^} (1) 

where Li denotes the contour length of IV The polymers can be either open or closed 
{y[ = r^. The topological states will be distinguished using the Gauss linking number: 

*^ r >> - h C C *<<«> * *'<«> ■ l £KrS (2) 

where xfiu Tj) takes integer values if Tj, Tj, i ^ j = 1, . . . , N, are closed trajectories. 
To describe the statistical mechanics of the polymers we define the configuration probability 
G{ m }({r}, {L}; {r'}, 0). This function measures the probability that the trajectories Tj fulfill 
the constraints: 

x(Ti, Tj) =mij i^j = l,...,N (3) 

and have extrema (in the open case) at the points and r^, % = 1, . . . , N. In our notations 
{m} denotes the n x n symmetric matrix of topological numbers with elements 

rriij = rriji for i ^ j ; ma = (4) 

while {r} = r 1; . . . , r^r, {L} = L%, . . . , L N etc. Starting from the path integral approach of 
Edwards 0, one obtains the following expression of G{ m }({r}, {L}; {r'}, 0): 
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G {m} ({r},{L}; {r'}, 0) 



+00 



N-l N 

nn 

i=l 3=2 

3>i 



2tt 



G { x } ({v},{L};{v'},0) (5) 



where 



G w ({r},{I};{r'},0)= Z>n(si) . . . / P riV (sjv)exp {- (Aq + A ev + A top )} (6) 



In the above equation the action 
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2af^Jo 

is that of a particle free random walk corresponding to a trajectory Tj, while 

a™ = ^ E / / ^K/ 3) ( r ^) - r i(4)) 



(7) 



(8) 



takes into account the excluded volume interactions. For convenience, v% has been defined 
as follows: 



v for i = j 
v%/2 for i 7^ j 
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(9) 



where the are coupling constants with the dimension of a volume. Finally, the topological 
constraints (|3]) are responsible of the topological term: 



N-l N 



4, = ^^^,,^ (10) 

i=l 1=2 

3>i 

coming from the well known Fourier transform of the Dirac 5— function: 



$( m i3 -X( r ^ r i)) 



2tt 



■exp {-i [rriij - xi^i, Fj)) Kj} 



(11) 



To express the configurational probability (|6]) in terms of fields, we introduce N gaussian 
scalar fields 0i(r), i = 1,...,N and N(N — 1) free C-S fields and B$n, with i = 
1, . . . , N — 1, j = 2, . . . , N and j > i. The action of the scalar fields is given by: 

,2 n 



(12) 



i,j=l 



while C-S one is given by: 

N-l N 



S CS = ^ /A E IX', • (V x Bg>) (13) 

^" J 1=1 3=2 

3>i 

The C-S theory will be quantized in the Landau gauge, where the fields A^l and are 
completely transverse. It turns out that the following relation is valid for the excluded 
volume term: 



- = Jvfa . ..V<f> N exp|-^ } - iJ2j ^(r^r) j (14) 



where 



Mr)= I 1 d Sl 5^(T-r t ( Sl )) (15) 
J o 

To treat the topological interaction it will be convenient to define the following linear com- 



binations of the fields and B^: 



c (1) = E« He) 



J=2 



N N-2 

j=3 j=l 



C W = E «f?)A$ + E i = 2, . . . , TV - 1 (17) 



and 



C (A ° = E /$)B« (18) 



i=i 

The parameters and depend on the C-S coupling constant k and on the matrix 
elements as follows: 

In terms of these fields C it is now possible to state the Gaussian identity: 

N ,L 



J VAVB exp^-iS CS -ipJ o * C«(r( Si ))dr( Si ) J = exp 



N-l N 

-i 



1=1 J=2 



(20) 
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Eq. (pOD, where we have used the compact notation: 



r N 

VAVB = / [] VA^VB^ (21) 

i<j=i 

relates the topological term (|K]) to the amplitude of N holonomies and will be central in the 
subsequent discussion. Substituting eqs. ( |2T)D and ([TJ]) in (Q), the configurational probability 
G {A}({r}, {L}; {r'}, 0) becomes: 

N 



G {A} ({r}, {I/}; {r'}, 0) = ([] Gfo L t] rj, 0)0,, C«)> w>{Am B> ( 22 ) 

i=i 

In the above formula the functions G(rj, L^; r^, 0\<f>i, C®) are given by: 

G(r i; L <; r^, O|0 i; C«) = j£* Pr^^exp |- ^ ds^r^)) - z ^ ds^s,) • C^(r^))} 

(23) 

where 

^X^)) = ^v%s l ) + i<p i (v i ) (24) 

Moreover, everything inside the bracket ( ){</>} {A} {B} mus t be averaged with respect 
to the fields {<p}, {A} and {B} by means of eqs. ( |T4j ) and (|20|). As we see from eqs. 
(p2f)-fl2~4D the polymer trajectories ri,...,rjv are completely decoupled before averaging 
over the auxiliary fields. Furthermore, each of the factors G{ri, L^; r' i: 0\(f)i, C^) appearing 
inside the bracket ( ) coincides formally with the evolution kernel of a particle 

subjected to random walk (in the Wiener sense of the path-integral formulation of statistical 
mechanics) and immersed in the electromagnetic field (ifa, C®). Thus G(i"j, Lf, r' { , O|0», C^) 
satisfies the pseudo-Schrodinger equation 0j: 

d a. 



<9L; 6 

The covariant derivatives Dj are given by: 



G(r h U- v>, 0|&, C«) = SiLiMn - rj) (25) 



D i = V + zC w i = l,...,iV (26) 
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At this point it is convenient to consider the Laplace transform of the configurational prob- 
ability G{ A }({r}, {L}; {r'},0) with respect to the polymer lengths [/J: 

r+oo r+oo ( J2_ 

G {x} ({r},{v'},{z}) = J q d Ll ...j Q rf ijv exp|-^^|G {A} ({r},{L};{r'},0) (27) 

Each variable Z{ plays the role of a Boltzmann-like factor governing the distribution length 
of the i— th polymer. Applying the above Laplace transformations to the expression of 
G {A}({r}, {L}; {r'}, 0) given by eq. (||), one finds: 



N 



G {x} ({v},{r'},{z}) = (II G '( r ^ r ^#^ c(i) )) Wl{ A M B} 



(28) 



i=l 



where the functions Gfa, r-; Zi\<fii, C^) are the Laplace transforms of Gfa, Lf, r' { , 0\4>i, C^] 
Therefore, they obey the stationary pseudo-Schroedinger equations: 



Zi ~ 6 T>i + ^ 



(29) 



The solution of eq. (|29| ) can be given in terms of second quantized fields fields ip*,ipi, 
z = l,...,JV: 



G(v i ,T f i ;z i \<j )i ,C®)=Z- 1 J 



(30) 



where F^i] represents the Ginzburg-Landau free energy of a superconductor in a fluctuating 
magnetic field: 



d 3 r 
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|D^| 2 + (zj + i^i)!^ 



(31) 



and Zi is the partition function of the system: 



(32) 



In order to eliminate the auxiliary fields 4>i, we exploit the method of replica ||20|| . Let 
ip*^, ipf* be a set of replica fields with oji = l,...,n, and i = 1,...,N, that form the 
n— pies = (if)}, . . . , ^r*) an d ^* — (V^* 1 ; • • • jV^ 71 *)- 111 terms of these fields the Green 
functions G(r i7 r-; C^) become: 



G(r 1) r' i ;^|0 i ,cW)= lim /wjCWrW e — 



(33) 



where the uJj's are arbitrary integers chosen in the range 1 < cJj < nj. According to the 
replica method, we also assume that the limit for rii going to zero commutes with the 
functional integrations over the C-S fields and the auxiliary fields fa. Substituting eq. ( |33| ) 



in (28) and performing the Gaussian integration over the fa, one obtains: 



N N 

hm /PAOTn^^ankrC^rW)] exp{-.4({*},{A},{B})} (34) 



i=l 3=1 

where the polymer free energy ^4({ 1 I'}, {A}, {B}) is given by: 

A({*},{A},{B}) = iS CS + 



N 

i=i 



|Dj\l/j| -|- 

o 



N „,0 



(35) 



«J=1 



From the above equation we see that the fields \&* have the dimension of a mass, which 
is not the canonical dimensionality of the scalar fields in three dimensions. To remedy, we 
introduce a new mass parameter M f2"I|| . It is nice to see that in this way it is possible to 
define an analog of the Planck constant: 

Ma 



(36) 



Performing the field rescaling 




(37) 



and working in units in which h — 1, the action (35) becomes: 



TV 



^({*'},{A},{B})=*5 C5 + E K r 



i=i 



+ 



* 2M 2 v°- 

^ 4J 'dW.| 2 |*' 12 



(38) 



where the covariant derivatives are defined as before by eq. ( p6|) and 

m\ = 2Mz, t (39) 

From now on, we will always use the rescaled action and the prime indices carried by 
the complex scalar fields will be dropped. 



III. THE TWO POLYMERS CASE 

In the two polymers case the matrix rriij of topological numbers is replaced by a single 
topological number m which describes the intersections of the trajectories r\ and T 2 . Its 
Fourier conjugate variable will be called A. It will also be convenient to put: 

Agj = AW Bgj = A( 2 ) (40) 

and 

a$ = 7i flg = 72 (41) 

In this new notation the action fl38D becomes: 

2 

A 2 (m, {A}) =i^J d 3 rA^ ■ (V x A< 2 >) + E jd 3 r [^*(-D 2 + m 2 )^ 



2 2M 2 ?; /• 

+ E — ^ / ^I^H^I 2 (42) 



a 2 

i,j=i a 
where 

D, = V + 2 7 ;A W (43) 

Accordingly, the configurational probability G{\}({r}, {r'}, {2;}) can be written as follows in 
the two polymers case: 

G A ({r},{r'},{,}) = 
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2 2 

lim /n [PA«P^P**] I] [#(r,)^r(rO] exp {-^(W, {A})} (44) 
Since we are interested in the effects of the topological interactions, we neglect the excluded 



volume interactions setting u£ = in the remaining of this Section. 



FIG. 1. Non-vanishing contribution to the one-loop effective potential due to the topological 
interactions. Solid lines represent the charged scalar bosons and dashed lines Chern-Simons vector 
fields. 

Tho study the above topological Ginzburg-Landau model at one loop order we exploit 
the method of effective potential |22| . To this purpose one has to compute only the diagrams 



of fig. ©. 



FIG. 2. Diagrams containing at least one vertex of the coupling A T ^ T ^/*. These contributions 
to the effective potential are identically zero in the Landau gauge. 

The other possible diagrams are in fact of the form given in fig. (^) and are suppressed 
in the Landau gauge ||22[| . First of all one needs to derive the Feynman rules starting from 



the action (f42"l). Here we just state the result: 

Air r>P 

(4 2 )(p)4 1 )(-p)) = »,v,p = 1,2,3 (45) 

K, p 

<C 1 (P)C J (-P)> = U = !> 2 (46) 
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(C 1 (Pi)^r(-P 2 )^(P 3 )) = 27^ (3) (Pi - P 2 + Ps)(P2), (47) 

(c i (Pi)^r(p 2 )4 i) (P3)4 i) (p 4 )) = t?m (8) (pi + p 2 + p 3 + p 4 ) ( 4 8) 

After some calculations, one obtains the following expression of the effective potential at 
one loop order: 

v.//(w) = i:™m 2 - ^i^n^r + ^ia^h^d 3 (49) 

In eq. (f49|) A is an ultraviolet cut-off with the dimension of a mass. Due to the fact that 
the microscopic scale of our model is given by the step length a, phenomena whose spatial 
extensions are smaller than a have no meaning [|7|, so that we can put: 

A ~ a' 1 (50) 

From eq. ( f4"9"D one sees that the total effect of the topological forces at one loop is to influence 
the strength of the repulsions among the polymers. In particular, topological forces do not 
affect the scaling behavior of the system. Indeed, only the last term appearing in the second 
member of eq. fl4~9|) may influence the critical behavior of the topological Ginzburg-Landau 
model fl38|) , but it is of higher order in the numbers of replicas n\ and n^- Therefore, it does 
not contribute in the configurational probability (|34]) , after the limit n^n-i — > is taken. 

The simplest observable that can be computed within the above approach is the second 
topological moment 

{m ] ~ rooG m ({r},M,{, } ) (51) 

Here G m ({r}, {r'}, {z}) is easily evaluated by taking the inverse Fourier transformation with 
respect to A of the configurational probability (f44|). This has been done in ref. ||17|| . In the 
case of the second topological moment one finds after some calculations: 

2 2 



(m 2 



}= lim /IJ [PA«P^P^*] J] [C(r^r(rD] exp{-A} 



i=l j=l 
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I dh^A^ ■ W 2 ) 2 + l -J rf 3 r|^ 2 | 2 A^ 2 > • A® 



(52) 



where 



A = i^-J d 3 rA (1) ■ (V x A (2) ) - J d 3 r^ 2 (D 2 - - J d 3 r^(A - m^ 2 (53) 



and 



= /n [x>A«2>* 4 2?*r] n [vf (r^rw; 



exp{-A} 



(54) 



Let us notice that in the action A$ the field A*- 2 - 1 plays the role of a Lagrange multiplier. As 
a consequence, after performing the integration in A^ 2 -* in eq. (0), the field A^ 1 -* becomes 
trivial, so that the path integral Z actually describes the fluctuations of two free polymers. 



¥i v; 



1 

2 



FIG. 3. This figure shows all the diagrams contributing to the second topological moment to- 
gether with the weights they appear in the expression of (m 2 ). 

The path integral appearing in the numerator of eq. (0) is instead more complicated. 
Nevertheless it is possible to show that, in the absence of excluded volume interactions, 
the only contributions to (m 2 ) come from the one loop diagrams of fig. (|3|). In this way 
the second topological moment can be exactly computed. Physically, one is interested to 
study the topological moments as a function of the length of the polymers. This requires a 
slight modification of the above approach to compute directly the configurational probability 
G{A}({ r }, {L}; {r'}, 0). The details of this calculation will be presented elsewhere. 
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IV. CONCLUSIONS 



In the previous Sections the statistical-mechanical problem of an arbitrary number of 
topologically entangled polymers has been mapped to a topological Ginzburg-Landau model, 
solving in this way a longstanding problem firstly pointed out in ref. |J. In our model the C-S 
fields acquire a physical significance, since they propagate the forces acting on the polymers 
due to the presence of the topological constraints. Of course, the topological forces should 
not distinguish a polymer from the other, a principle that is apparently violated by our choice 
of coupling constants in eq. (EH^). However, explicit calculations show that this symmetry is 
not spoiled and actually puts severe constraints on the number of possible Feynman diagrams 
to be computed in the perturbative approach. For instance, one recognizes from eq. fl49[) 
and fig. (0) that the contributions to the one loop effective potential and to the second 
topological moment are completely symmetric in the fields and ^2- A related 

problem is the appearance of a spurious parameter in the polymer free energy fl38|), namely 
the C-S coupling constant k. Also in this case it is easy to see that the difficulty is only 
apparent, since it is possible to remove the dependency on k with a simple rescaling of the 
fields Ag. 

In our construction the topological forces acting on each couple of trajectories Tj and 
Tj are propagated only by the two fields A^l and B^l. For that reason, the calculations 
performed in the case of two polymers can be extended also to an arbitrary number of 
polymers in a straightforward way. In this way it is possible to conclude from the analysis 
of Section III that the topological interaction does not change the critical behavior of the 
system at least in the one loop approximation. 

As already mentioned in the Introduction, the Gauss linking number is a rather poor 
knot invariant. However, the field theoretical formulation of topologically linked polymers 
established here shows a natural way to include higher order topological interactions. This 
is achieved by adding to the action (|3^) a new contribution of the form: 
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N-l N . 

A int _ \ - \- / rfir^ up ii ijk A (i) A (i) A {k) +fi ijk A {i) A (i) R (A:) 4- "l (55) 

m)y (n)p Imn (m)v (n)p 1 ' " 'J \ J 

i,j,k=l i,m,n = 2 

i<Z;j<m;fc<n 

Each term appearing in the above sum is topological and does not vanish apart from special 
combinations of the indices k, I, m, n, in which two or three C-S fields coincide. Clearly, 
interactions as those given in eq. fl55[) generate corrections to the configurational probability 
G{\}({r}, {L}; {r'}, 0) that contain higher order link invariants. In perturbation theory 
these corrections can be computed order by order with techniques similar to those already 
exploited in the case of nonabelian C-S field theories [0 , jl9) . To enforce the new topological 
constraints through Dirac 5— functions as in eq. ([|), one has to introduce additional Fourier 
variables besides the parameters Ajj's. The possibility of defining the coupling constants 
limm ^imn e ^ c - m terms of these variables in a way that is suitable to impose relations on 
higher order link invariants is still an open problem and deserves more investigations. 
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